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ABSTRACT 



In this thesis a clutter voltage model for scattering from the sea surface is developed. 
A model for the scattering from a whitecap and a wave breaking occurrence model t^re 
combined to simulate the back scattered signal from one radar resolution cell. The 
simulation performed obtained the probability density function of sea clutter under dif- 
ferent assumptions of wind velocities and wave breaking conditions. This model incor- 
porates some measured quantities such as the mean clutter voltage and the correlation 
time as parameters. The probability density function depends on the parameters of this 
model. The obtained probability density functions do not confirm to any familiar simple 
density function. 
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I. INTRODUCTION 



In a radar system, the returned signal from a target is contaminated with clutter. 
The decision of a target presence is made based on the statistical description of targets. 
The clutter voltage in the envelope detector output is considered in this thesis and is 
treated as a random variable. Note that only the magnitude information but not phase 
is preserved by envelope detection. Therefore, in this thesis, "clutter voltage" will always 
mean the magnitude only. 

The receiver output power is proportional to clutter cross section. The clutter volt- 
age is proportional to the square root of clutter cross section. 

In a radar receiver, components unavoidably generate thermal noise. The detector 
voltage output is contaminated with thermal noise. Thermal noise is usually considered 
as following the Gaussian distribution. 

Figure 1 demonstrates the principle of threshold detection which is widely used in 
radar systems. Curve P 0 (y) is the probability density function of the voltage output when 
there is no target present. The mean of this probability density function is zero. Curve 
P,(y) is the voltage when there is a target present ( Reference 1). 

The decision threshold is determined by Neyman-Pearson criterion, which maxi- 
mizes detection probability while keeping a constant false alarm rate. There are two 
kinds of errors. One is the false alarm. The false alarm probability is the area above the 
threshold under curve P a (y). The other is the probability of missed detection. Its prob- 
ability is the area below the threshold under curve P,(y). For a constant false alarm rate 
detection using the Neyman-Pearson criterion, the density function without the presence 
of a target determines the threshold level. The density function with the presence of 
targets determines the missed-detection probability. The same applies in the detection 
of a target in clutter. Clutter is the unwanted back scattered signal from the sea surface. 
It is a coherent RF signal similar to the target return. Thus the integration techniques 
which are designed to increase the signal to thermal noise ratio can not be used against 
clutter. The clutter power is proportional to the illuminated area of a resolution cell of 
the radar. 

Figure 2 is the geometry of a resolution cell. The illuminated area is the product of 
range, beam width, range cell and the secant of the grazing angle. The received power 
is directly proportional to the clutter cross section, sigma. The proportionality constant 
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Pi(y) 




Figure 1. Decision Diagram in Gaussian Noise. 

is a function of the space attenuation and antenna gain. Conventionally, the clutter 
reflectivity, which is the clutter cross section per unit area, is used. Note that the 
reflectivity, being a ratio of the cross sectional area to the illuminated area of the sea 
surface, is a dimensionless constant. At X-band if the range cell is greater than 20 in 
(Reference 1), the radar is considered as lo\V resolution radar. 

The radar signal is usually narrow band. Any narrow band signal may be expressed 
as x(t)coscot + y(t)sincot, where x(t) and y(t) are the 1 channel and Q channel signals. For 
a low resolution radar, the signals are scattered from a fluctuating surface of large illu- 
minated area. By the central limit theorem, both the I channel and Q channel signals are 
Gaussian. The envelope detector output voltage, then, is Rayleigh distributed and is 
called Rayleigh clutter. 
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Figure 2. Geometry of Resolution Cell. 

Modern radars have higher resolution. The reason is that by decreasing the illumi- 
nated area, they increase the signal to clutter ratio. The power scattered from a target 
remains constant if the target is in a resolution cell. The power scattered from the sea 
surface is reduced since the resolution cell is smaller. 

It is observed though that, while the total back-scattered power from the sea surface 
is reduced, the peak power is not decreased and the clutter signal becomes spiky in ap- 
pearance. Thus the false alarm rate is not reduced even though the clutter power is 
lower. Observations associate such spikey nature with the presence of whitecaps (Ref. 
1 ). 

Measured data show the reflectivity of breaking water to be about 1 to 1.5 while that 
of unbroken water is about -40 dB below unity (Ref. 1). Hence, when breaking waves 
are present, their back scattering is the major components of the sea clutter. 
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Furthermore, since the number of breaking waves within a resolution cell is limited, 
the central limit theorem does not apply to the clutter of a modern, high resolution ra- 
dar. This may be the reason for the non-Rayleigh distribution observed in many meas- 
urements. 

Figure 3 is one example of the non-Rayleigh sea clutter. Data obtained by NRL 
show that the distribution has a longer tail compared to the Rayleigh distribution (ref- 
erence 8), which means there is a higher probability of stronger sea clutter. The simu- 
lation of this phenomenon has been carried out by Lewis and Olin (Reference 1) and 
will be discussed in chapter 2. 

In the investigation of clutter distribution, many effects should also be considered. 
The reason that sea return is random is that the surface is fluctuating. Many researchers 
describe the water surface using the slope and height of the waves as the random vari- 
ables. In terms of these random variables, the conditions for wave breaking can be stated 
as some breaking criteria (Reference 2) derived from ocean wave dynamics. 

Chapter 3 provides a model for the back scattering from the sea surface. Since 
scattering is very different between broken and unbroken water waves, a breaking oc- 
currence model is developed. Back scattering from unbroken water and from breaking 
waves are simulated separately and then combined afterward according to the positions 
and sizes of white caps at the sampling time. 

Chapter 4 is the simulation of total back- scattering signal from the sea. In each 
simulation run of different parameters, 20,000 independent samples are calculated. They 
are used later to find the PDF of clutter voltage. 

In Chapter 5, the PDF of the clutter voltage is calculated after combining back- 
scattering from whitecaps and unbroken water. The clutter voltage is truncated and di- 
vided into 500 equal intervals. The number of samples located in each internal divided 
by 20,000 is the probability of the clutter voltage occuring in this interval. 

The simulated PDF is obtained by dividing the clutter voltage probability with the 
width of the interval. 

The simulation result shows that clutter voltage from unbroken water fits well with 
a Rayleigh distribution if the randomness factor to be defined in Chapter 2 is 1. Other- 
wise it may be fitted with the Weibull distribution. This compares favorably with some 
measured data (Reference 3). 

Clutter distribution from broken water has a new form. This thesis suggests that its 
shape depends on the sea state and the randomness factor of the white caps. 
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Figure 3. Experimental sea clutter cross section data taken by NRL on X band 
radar, 0.02 microsecond pulse, vertical polarization and 4.7 degree grazing 
angle, and Rayleigh distribution for comparison. 
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II. SCATTERING FROM A WHITECAP OR AN UNBROKEN SURFACE 



A. BASIC THEORY 

Calculating the backscattered signal from the sea surface involves a better under- 
standing of ocean dynamics, beyond that which is currently available. 

The goal for this thesis is to find a framework upon which measured parameters can 
be incorporated to provide a statistical description of the sea clutter. A model developed 
by Olin and Lewis is adopted to calculate the scattering from growing whitecaps. This 
model and its implemention are discussed in this chapter. In later chapters, this simple 
model is extended to include the spatial and temporal statistics of the breaking waves to 
obtain the PDF of the sea clutter. With the extended model, the spatial and temporal 
correlation of the clutter voltage can be investigated later. 

Lewis and Olin's model is for the back scattered signal from a wave crest, including 
the white cap. In this model, a sea wave is divided into identical zones of a quarter radar 
wavelength each. Within each zone, the scattered field from the sea surface is roughly 
in phase and is considered to be due to a source located at the center of this zone. The 
relative magnitude of the back scattered signal from this zone is 1-kX and the phase is 
2;rkY. Both X and Y are random variables uniformly distributed over 0 and 1. Therefore 
when k equals 1, the magnitude of the back scattering can van' over the complete range 
of 0 to 1 while the phase of the signal can vary over the complete range of 2n. This 
parameter k is called the ‘'randomness factor'*and the back scattering is completely ran- 
dom when k = 1. When k does not equal to unity, a nonfiuctuating component in the 
back scattered signal is introduced. Thus k can be considered as a measure of the 
S'oughness'^ or lack of) of the wave. 

Assume a radar is operating in X band with a 3 cm wavelength. Each quarter 
wavelength zone is .75 cm. The growth rate of a w'hite cap is observed to be 1 m/sec. 
The time for each white cap to grow one zone larger is 0.75 milliseconds. This approx- 
imately agrees with the 1 millisecond decorrelation time observed. The w'hitecap grow's 
linearly from length zero to 0.96 meter, when the number of zones grow's from 0 to 128. 

Consider only the case w'hen the sea waves are moving toward the radar. The back 
scattered signal referred to the back edge of the w'hite cap is 
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e m = s rn ex P -j 4 y C(W- 1) y +4J 



( 1 ) 



where : 

5 m < 1 = equivalent zonal reflection coefficient. 

0 <d m <,-^r = location of effective scatterer within the zone. i 

S m = ( \-kX m ) (2) 

^ = T + 4 2 0 (3) 

where : 

0 < X m < 1 = random component in amplitude. 

0 < Y m < 1 = random component in phase. 

0 < k < 1 = fraction of random component. 

By using these expressions and shifting the initial reference surface by the scat- 
tered signal can be written: 

= (-ir’d - kX m ) exp —jnk(Y m — ± ) (4) 

IV 

If the range extent of the whitecap is W, then the number of zones is N = 4 —y , 
and the total backscattered signal f A - is given by 



a- 

r . v : = Yji- 1 T ' 1 ( !-**«) exp -jnk( Y m -\) (5) 

m= 1 

where k is the randomness factor of the back scattered signal. It is assumed to be con- 
stant from zone to zone. Note that only the relative magnitude of the scattered signal 
is considered. The mean clutter voltage is determined from experimental measurement. 

In order to model the growth of a white cap, D is introduced as the internal decor- 
relation time. Hence in each zone, the phase and amplitude of the scattered signal will 
be another pair of random numbers after the white cap grows by D quarter wavelength. 
For example, if D is 5, then in this thesis the phase and amplitude of scattering from 
effective scatterers in every zone of a growing whitecap will be another pair of random 
variables after 5 times 0.75 millisecond. 
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In the computer simulation of scattered signals from breaking waves, phases and 
amplitudes are generated from a uniform random number generator. The growth of a 
whitecap is described by increasing its zone number. 

B. COMPUTER SIMULATION 

The algorithm to simulate back scattering from one growing breaking wave is : 

1 set whitecap size from 1 to 128, 

2 call the uniform random number generator to 

get 2 random numbers for the phase and the amplitude for 
each zone, 

3 sum the back scattered signal from all zones at 
the sampling times, 

4 repeat for each sampling time, and then 

get the maximum value of samples and normalize 
all samples. 

Since the scattered signal is time-dependent, it is not a stationary' process. A 
FORTRAN subroutine is listed in Appendix A to simulate the scattered signals for one 
breaking wave growing process. The output of this subroutine is an array containing 128 
scattered samples. The mean of each set of samples is different. But according to Lewis 
and Olin's measurement, the reflectivity of a breaking wave is about 1 to 1.5 
(dimensionless or M**2/M**2). We choose a value of 1.5 as the reflectivity of the 
whitecaps. The mean clutter cross section is known to be the reflectivity times the 
whitecap area. 

The clutter voltage is 1.224745. So in the simulation, the 128 samples are first nor- 
malized by the maximum sample value. Then they are scaled by 1.224745 as the mean 
of the samples. Figure 4 is a set of data measured by Lewis and Olin. It shows the 
clutter cross section from a single resolution cell. A modulation imposed on the signal 
is evident which is attributed to the intermittancy of the breaking waves (Reference 3). 

Figure 5 is a realization of this model with k= 1 and D= I from one breaking proc- 
ess. Figure 6 is a realization with D= 1 and k = 0.8, Figure 7 with D= 1 and k = 0.1. 
there is little modulation on samples of small zone numbers. Figure 8 is a simulation 
with D= 1, k=.06 and modulation becomes more obvious. Figure 9 is with k= 1 and 
d = 5 and shows wider fluctuation. 
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Figure 4. Time history of clutter cross section measured by Olin and Leu is 



WHITE CAP SIMULATION 



Figure 5. 




Simulated scattered signal from uhitecap, D = 1, k = 1. 
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Figure 6. 



WHITE CRP SIMULATION 




Simulated scattered signal from uliitccap, D = 1, k = 0.8 
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WHITE CRP SIMULATION 



Figure 7. 




Simulated scattered signal from uhitecap, D = 1, k = 0.1. 
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WHITE CRP SIMULATION 




Figure 8. Simulated scattered signal from nliitecap, D = 1, k = 0.06. 
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WHITE CAP SIMULATION 




Figure 9. Simulated scattered signal from whitecap, D = 5, k = 1. 

The same model is also used in simulating the scattering from the unbroken water. 
The dilference is that the area of the unbroken water is the whole resolution cell. This 
is a rectangle of 225 square meters. The range cell is 15m for a pulse width of 100 
nanoseconds. The range cell contains 2000 zones, each zone has a length of 0.75 cm. 

i 

Each clutter sample is the sum of scattered signals of 2000 terms. In the simulation, the 
uniform random number generator is called 2000 times to generate random phases and 
amplitudes. There is no growing effect in the unbroken water. 



C. AUTOCORRELATION 

Because of the growth of whitecaps, the clutter voltage is not a stationary process. 
The number of terms is changing, so it is not expected to have the same autocorrelation 
function. 
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Figures 10 to 16 are autocorrelation functions with different randomness factors and 
internal decorrelation times. 
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HH1TE CHP RUTOCORRELflTION 
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Figure 12. Autocorrelation function of scattering from uhitecap, D = 1, k = 0.5. 
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WHITE CAP AUTOCORRELATION 




Figure 14. Autocorrelation function of scattering from nhitecap, D = 1, k = 0.2 
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WHITE CAP AUTOCORRELATION 




Figure 15. Autocorrelation function of scattering from uliitecap, D = l,k = 1 
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WHITE CAP AUTOCORRELATION 




Figure 16. Autocorrelation function of scattering from whitecap, D = l,k = 1 , 
seeds of random number generators are different from seeds in Figure 15 



D. PROBABILITY DENSITY FUNCTION 

It is well known that the clutter from a large observed area obeys the Rayleigh dis- 
tribution. The scattering model developed -above is used to simulate the scattering from 
unbroken water. It is clear from the results of this section that the probability density 
function of such clutter cross section is an exponential distribution. Thus the model is 
verified for this particular case. 

There are 2000 independent samples simulated to calculate the probability density 
function of clutter cross section from an unbroken water. The result shows that when 
the randomness factor k is 1, the probability density function of the clutter cross section 
is an exponential distribution, while the clutter voltage is a Rayleigh distribution. 
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Figure 17 and 18 arc PDFs with k= 1. The dots are calculated PDF from the sam- 
ples. The average of these samples is used to plot the theoretical distribution of clutter 
cross section, which is the solid line. They arc both exponential distributions. So the 
scattering model developed by Olin and Lewis is valid for unbroken water if the ran- 
domness factor is 1. Figures 19 and 20 are RCS distributions with a randomness factor 
of 0.3. 



RCS WITHOUT WAVE BREAKING 




Figure 17. Probability density function of sea clutter cross section, randomness 
factor = 1. 
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RCS WITHOUT WAVE BREAKING 

s 




normalized rcs 



Figure 18. Probability density function of sea clutter cross section, randomness 
factor = 1. Seeds of random number generators are different from seeds 
in Figure 17. 
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RCS WITHOUT WAVE BREAKING 




Figure 19. Probability density function of sea clutter cross section, randomness 
factor = 0.3. 
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RCS WITHOUT WAVE BREAKING 



i 




Figure 20. Probability density function of sea clutter cross section, randomness 
factor = 0.3. 

Figures 21 and 22 arc PDFs with randomness factors of 0.01. They deviate from 
an exponential distribution more significantly. Because the scattering model for unbro- 
ken water is verified to agree with the scattering of Rayleigh clutter if the randomness 
factor is 1, it is not necessary to sum up all 2000 terms for every sample. The Rayleigh 
random generator can be used and its mean is to be scaled by the unbroken water 
reflectivity. 
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RCS WITHOUT WAVE BREAKING 




NORMAL I ZEO RCS 



Figure 21. Probability density function of sea clutter cross section, randomness 
factor = 0.0 1. 



26 



RCS WITHOUT WAVE BREAKING 

o • 




Figure 22. Probability density function of sea clutter cross section, randomness 
factor = 0.01. 
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III. A SEA SURFACE MODEL 



In Chapter 2, the scattering model simulates back scattered signals from the unbro- 
ken water surface and the individual breaking wave separately. Since, at any sampling 
time, there will be some breaking waves in the radar illuminated area, the locations and 
sizes of whitecaps must be modeled and combined with the unbroken sea surface to 
calculate the total backscattered signal. 

In this chapter, a sea surface model is developed to simulate the occurrence of 
whitecaps in the resolution cell. 



A. OCCURRENCE OF WHITECAPS 

It is known that the occurrence of a white cap is due to the interaction of wind and 
long ocean waves. When the wind blows, the capillary waves are formed and superim- 
pose themselves on the long, faster waves. When the combined velocity of water waves 
on top of the long waves attains a kinetic energy beyond the surface tension potential, 
the water breaks away from the surface. This process is analogous to shot noise; when 
a kinetic energy threshold is reached, an electron breaks away from the metal surface and 
creates shot noise in the circuit. 

Thus to describe the temporal occurrence of a breaking wave at a location, the 
Poisson distribution is adopted. Since the occurrence and kinetic energy of the capillary 
wave depend on the wind velocity, the average rate of the occurrence of the wind ve- 
locity will be given later in this chapter. The position of each breaking wave is assumed 
to be two dimensional, uniformly distributed, over the resolution cell. 



B. COMPUTER SIMULATION 

In this section, a FORTRAN subroutine listed in Appendix B is discussed The pur- 
pose of this subroutine is to simulate the occurring times and locations of breaking 
waves. The total observing time interval was 750 seconds. A Poisson random number 
generator is used to simulate a sequence of breaking wave occurrence times. 

A uniform random number generator was used to simulate the locations of breaking 
waves. The illuminated area was 225 square meters and was divided into nine square 
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sub-areas. In each sub-area the occurring times and locations of whitecaps were simu- 
lated independently. 

The procedure of this simulation is to call the Poisson random generator to get a 
sequence of starting times over the entire observation period. Each starting time is the 
beginning of one breaking wave. For each whitecap, the uniform random number gen- 
erator is called twice to get the X position and the Y position for each whitecap. This 
procedure is done nine times independently for the nine subareas. The nine sequences 
of occurring times are sorted and put into an array. Two other arrays are used to store 
the X and Y positions of the whitecaps corresponding to each of the occurring time se- 
quences. 

At any sampling time, the existing whitecaps are those with occurring times leading 
the sampling time by no more than 0.96 seconds. The time differences are the sizes of 
existing whitecaps since the growth of whitecaps is assumed to be linear. 

Figure 23 to Figure 34 are consecutive samples of the observed area over one hun- 
dred samples, with a randomness factor of 0.9. The sampling time increment is 0.2 sec- 
onds. The time period is from t = 25.4 seconds to t = 27.6 seconds. 
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Figure 23. Whitecaps occurrence realization, total observation time 100 seconds, 
randomness factor = 0.9, sample time at 25.4 second. 
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Figure 24. Whitecaps occurrence realization, total observation time 100 seconds, 
randomness factor = 0.9, sample time at 25.6 second. 
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Figure 25. Whitecaps occurrence realization, total observation time 100 seconds, 
randomness factor = 0.9, sample time at 25.8 second. 
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Figure 26. Whitecaps occurrence realization, total observation time 100 seconds, 
randomness factor = 0.9, sample time at 26.0 second. 



33 



WHITE CAP MODEL Klo 

26.200 





0.8000 




0.7708 

0.2556 




0.1633 




0.5696 


0.2122 



0.0 5.0 10.0 15,0 
X POSITION im* 



Figure 27. Whitecaps occurrence realization, total observation time 100 seconds, 
randomness factor = 0.9, sample time at 26.2 second. 
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Figure 28. YVhitecaps occurrence realization, total observation time 100 seconds, 
randomness factor = 0.9, sample time at 26.-1 second. 
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Figure 29. Wliitecaps occurrence realization, total observation time 100 seconds, 
randomness factor = 0.9, sample time at 26.6 second. 
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Figure 30. Whitecaps occurrence realization, total observation time 100 seconds, 
randomness factor = 0.9, sample time at 26.8 second. 
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Figure 31 
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Whitecaps occurrence realization, total observation time 100 seconds, 
randomness factor = 0.9, sample time at 27.0 second. 
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Figure 32. Whitecaps occurrence realization, total observation time 100 seconds, 
randomness factor = 0.9, sample time at 27.2 second. 
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Figure 33. Whitecaps occurrence realization, total observation time 100 seconds, 
randomness factor = 0.9, sample time at 27.4 second. 



40 



WHITE CAP MODEL So 

27.600 



2 

o 



to 

a 

Q_ 



0 . 26 - 4 1 




0.9509 


0.6211 




0.5900 


0.2270 






0.4492 


0.7004 


0.0000 



S.O 10. 0 

X POSITION (Ml 



15.0 



Figure 34. Whitecaps occurrence realization, total observation time 100 seconds, 
randomness factor = 0.9, sample time at 27.6 second. 



Note that all whitccaps remain at the same locations with their sizes increasing lin- 
early with time. When the size of a white cap reaches its full size of 12S zones, it disap- 
pears. 

< 

The rate parameter of the Poisson process are related to the sea state. Monahan 
and O'Muircheartaigh (Reference 5) conclude empirically that the whitecap coverage are 
related to the wind speed by the following equations: 

)V = 2.95 x 10~ 6 x 17 3 ’ 52 (6) 

or 

\V = 3.84 x 10~ 6 x U 3Al (7) 
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where U is wind speed in meters per second, and W is the fraction of the sea surface 
covered by white caps. 

The first equation is derived by the least-squares fitting of measured data. The sec- 
ond equation was derived using the robust biweighted fitting from the same data set. 
There are many other mechanisms afFecting the local wave energy density, such as wave 
current interactions at the strong current boundaries, sea surface boundary stability and 
the pre-existing sea state (Reference 6). 

At every sampling time, all sizes of the existing whitecaps are known. So the aver- 
age coverage percentage of the sea surface can be computed over all samples. 

The motion of wave is not modeled, because the velocity of the wave is assumed to 
be constant. Over any time delay, the difference in locations of existing whitecaps is 
constant. Then the motion contributes to an overall phase shift which does not affect 
the total scattered signal. 
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IV. TOTAL BACK SCATTERING SIGNAL 



In this chapter, the complete simulation program is developed to calculated a suffi- 
cient number of back scattered signals from the illuminated area containing whitecaps. 
Then the probability density function of clutter voltage is calculated as a function of 
various parameters. 

In previous chapters, a whitecap occurrence model was used to simulate the wave 
breaking process in one resolution cell. The illuminated area is 15 m by 15 m. The entire 
observation time is 750 seconds. A Poisson random number generator generates a se- 
quence of starting times and a uniform random number generator generates the lo- 
cations of each whitecap. These locations are the reference locations for the 
backscattering simulation and are assumed to be static since the wave motion contrib- 
utes no relative phase shift in the scattered signal, if the velocity of the wave is constant. 
It is assumed that even’ whitecap has a full size of 0.96 m, a width of 1 m, and is 
breaking toward the radar receiver. 

The size of a whitecap is assumed to be growing linearly at a rate of 1 m/sec. To 
calculate the total scattered signal at any time, the back scattered field from the unbro- 
ken water surface has to be combined with the contributions from the white caps over 
the illuminated area. 

These white caps are at different stages of breaking and have different sizes. Each 
size determines the sample number in the corresponding 128 scattering samples over the 
life of a white cap. The reference location of the unbroken surface is chosen to be the 
outer boundary of the resolution cell. 

The reference location of the scattered field from a breaking wave is where the wave 
breaking begins. So in summing these scattered fields in complex numbers, it is neces- 
sary to align them to the same outer boundary of the resolution cell. This is done by 
shifting the phase of each scattered signal according to the distance of the whitecap to 
the outer boundary of the resolution cell. 

Before summing the scattered signals, since they are all normalized sequences of 
samples, the mean strength of the unbroken water is scaled by 0.15 m. The reason is 
that according to actual observations, reflectivity of the unbroken water is -40 db below 
unity, or 0.0001. The clutter cross section is then 0.0001 times the area illuminated. So 
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the RCS due to the unbroken surface is 0.0225 m**2. Taking the square root, the clutter 
voltage is a normalized 0.15. 

The scattered signal from whitecaps is normalized at 1.224745 because a breaking 
wave's reflectivity is about 1.5. The area of one whitecap is 1 square meter, so each one 
has a clutter cross section 1.5 square meter. The clutter voltage is the square root of 
clutter cross section, so it is normalized at 1.224745 (Reference 1). 

After shifting the phase, all scattered signals are summed to give the total back 
scattered signal. Its amplitude is the clutter voltage. The clutter cross section due to the 
particular normalization adopted in this thesis, is the square of this voltage, in square 
meters. In every run of this simulation, 20,000 samples are calculated. 

A voltage range is selected that can cover all the samples. Dividing the range into 
500 intervals and counting the number of samples locate in each interval and then di- 
viding by 20,000, one gets a result which is the probability of the occurrence of clutter 
voltage within this interval. 

The probability density at the center of each interval is the probability of the ran- 
dom variable within that interval divided by the width of the interval, if the interval is 
small enough. 



44 



V. PROBABILITY DENSITY FUNCTION 



In this chapter, probability density functions of different randonuiess factors and 
internal decorrelation times are discussed. 

The parameters in the simulation developed above include an internal decorrelation 
time in the growing process of whitecaps, a randomness factor representing roughness 
of whitecaps, and a Poisson process rate parameter, which determines the number of 
breaking waves. When wind speed is very small, a small occurrence rate causes no wave 
breaking over the entire observation time. 



PROBABILITY DENSITY FUNCTION 
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Figure 35. Comparison of PDF between simulation and Rayleigh distribution. rate 
= 0.0001, observing time = 750 seconds, total whitecaps = 0. 
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Figure 35 shows the comparison between simulated probability density function of 
tire clutter voltage and theoretical Rayleigh distribution when the water is unbroken. 
The dotted line is result of the simulation; the solid line is the Rayleigh distribution 



with the mean equal to the simulated result. 
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Figure 36. Simulated PDF of clutter voltage, observing time 750 seconds, total 
nhitecaps is 74, Randomness factor = 1. 



When the wind speed increases, the rate parameter of the Poisson process increases. 
Figure 36 is the result of 74 breaking waves occurring over the entire observation time. 
The mean clutter voltage is 0.5692 , while the mean clutter voltage in Figure 35 is 0.1501. 
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Figure 37. Simulated PDF of clutter voltage, observing time 750 seconds, total 
whitecaps is 74, Randomness factor = 0.25. 



This increase in the mean clutter voltage is due to contributions by the breaking 
waves. This can be observed from the higher tail in the probability of Figure 36, com- 
pared to that of Figure 35. Notice the total probability in this simulation up to a clutter 
voltage of 6 is 0.9998. This means there iS only a negligible number of samples greater 
than 6. Note that the maximum sample value included is 6.1831. 
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Figure 38. Simulated PDF of clutter voltage, observing time 750 seconds, total 
whitecaps is 78, Randomness factor = 1. 



Figure 37 uses a simulation with the same seeds for random generators as those used 
in the simulation of Figure 36. In both cases there are 74 breaking waves occurring over 
the entire observation time. The only difference is that in Figure 37 the randomness 
factor of the whitecaps is 0.25 while the rahdomness factor in Figure 36 is 1. The con- 
tributions to the probability density function from the unbroken water and the whitecaps 
are similar. But in the conjunction part of the probability density function of Figure 37, 
the 0.25 roughness factor causes a sharper notch. 
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PROBABILITY DENSITY PUNCTION 




Figure 39. Simulated PDT of clutter voltage, observing time 750 seconds, total 
uliitecaps is 78, Randomness factor = 0.25. 



Figure 38 is a different run of the simulation, with 78 whitccaps occurring. The 
roughness factor is 1. Figure 39 uses the same seeds for random generators as those used 
in Figure 38, with the randomness factor of the whitccaps taken as 0.25. The difference 
Figure 38 and Figure 37 is similar to that between Figure 36 and Figure 37. 
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Figure 40. Simulated PDF of clutter voltage, observing time 750 seconds, total 
uhitecaps is 696, Randomness factor = 1. 
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PROBABILITY DENSITY FUNCTION 




Figure 41. Simulated PDF of clutter voltage, observing time 200 seconds, total 
uhitecaps is 1452, Randomness factor = 1. Solid line is Rayleigh dis- 
tribution with the simulated sample mean. 



In Figure 40. 696 whitecaps appear over the simulation period. The peak corre- 
sponding to the Rayleigh distribution of Figure 35 diminishes further because there arc 
stronger contributions from the breaking wives. The tail of the clutter PDF in this case 
becomes more pronounced. The maximum clutter voltage sampled increases to 7.831 , 
The mean clutter voltage increases to 0.9812. The total probability up to the clutter 
voltage of 4 in this figure is 0.985. The peak probability density at zero clutter voltage 
to 2.1 . 
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PROBABILITY DENSITY FUNCTION 




Figure 42. Simulated PDF of clutter voltage, 1000 scattering sample from whitecaps 
only, Randomness factor = 0.1. 



Figure 41 is a run of simulation with a greater rate parameter. The Poisson process 
generates many more breaking waves over the entire observation period. The observa- 
tion time is 200 seconds. There are 1452 whitecaps occurring. The peak of the original 

< 

Rayleigh distribution of Figure 35 caused by the scattering from the unbroken water 
now vanishes completely. The entire probability density function is dominated by the 
scattering from whitecaps. It is similar to the Rayleigh distribution with a mean of 
3.7286. 

The solid line in Figure 41 is the Rayleigh distribution with the simulated sample 
average as its mean. 

In Figures 36 and 38 the probability density functions show a sharper notch in the 
conjunction between contributions from the whitecaps and the unbroken water com- 
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pared to Figure 37 and 39. The differences are due to the fact that the randomness 
factor of whitccaps for Figure 37 and 39 is a smaller number. The probability distrib- 
ution of the scattering from growing whitccaps should be further investigated. Thus, a 
program is designed to collect samples from the breaking waves only. 



PROBABILITY DENSITY FUNCTION 




Figure 43. Simulated PDF of clutter voltage, 1000 scattering sample from whitecaps 
only, Randomness factor = 0.08. 

« 

Figure 39 is the probability density function of normalized clutter voltage calculated 
from 10.000 independent samples from whitecap. The randomness factor is 0.1. The 
reason for using normalized clutter voltage is that it is necessary to investigate the shape 
of probability density function when the randomness factor is small. In this figure, there 
is a ripple in the PDF. If the contribution in PDF from unbroken water overlaps the 
normalized clutter voltage from 0 to 0.3, then depending on the phase of the signals, the 
ripple at 0.3 may cause a sharp notch in the probability density function of the combined 
signals. 



53 



PROBABILITY DENSITY FUNCTION 
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Figure 44. Simulated PDF of clutter voltage from whitecaps only, Randomness 
factor = 0.06. 
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PROBABILITY DENSITY FUNCTION 




NORMALIZED RCS 



Figure 45. Simulated PDF of clutter voltage from whitecaps only, Randomness 
factor = 0.04. 
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Figure 47. Simulated PDF or clutter voltage from uliitecaps only, Randomness 
factor = 0.005. 
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PROBABILITY DENSITY FUNCTION 




Figure 48. Simulated PDF of RCS from white caps only, k = 0.6. 

Figures 42 to 47 are probability density functions of the clutter voltage scattered 
from growing whitecaps, with the randomness factors 0.1, 0.08, 0.06, 0.04, 0.02, and 
0.005 respectively. Figure 48 is the simulated PDF of RCS with k = 0.6. Note that the 
PDF's have two peaks and their peaks becomes farther separated when k decreases. 
Meanwhile the peak values of the probability density functions increase while the ran- 
domness factor decreases. For a randomness factor of 0.005, the probability density 
function is zero between 0.15 and 0.9 normalized clutter voltage. 

It can be concluded that the sharp notches in Figures 37 and 39 are caused by a 
small randomness factor. 

In Figures 49 and 50, the same set of seeds of the random number generator is used 
for the simulation, with randomness factors of 0.1 and 0.05 respectively. 
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Figure 49. 




Simulated PDF of clutter voltage, observing time 750 seconds, total 
uhitecaps is 78, Randomness factor = 0.05. 
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Figure 50. 
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Simulated PDF of clutter voltage, observing time 750 seconds, total 
vvliitecaps is 78, Randomness factor = 0.1. 
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Figure 51. 



Simulated PDF of clutter voltage from uhitecaps only, Randomness 
factor = 0.005, internal decorrelation time = 2. 
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Figure 52. Simulated PDF of clutter voltage from nhitecaps only, Randomness 
factor = 0.005, internal decorrelation time = 5. 

Figure 51 shows the results of simulation with an internal decorrelation time of 2. 
Figure 52 is the results of simulation with an internal decorrelation time of 5, with the 
same seeds of random number generators. It shows little difference when the internal 
decorrelation time is changed by this amouht. 
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VI. SUMMARY 



Many measurements have been made to analyze the probability distribution of sea 
clutter. Several probability density functions have been used by other researchers to fit 
those measured data. To date, none of them appears satisfactory'. A different approach 
is proposed to construct a model with a few physical parameters which describe the sea 
surface and then to obtain the clutter voltage in terms of simulation. Those physical 
parameters will be determined by fitting the simulated clutter voltage with measured 
data. 

In this thesis, the probability density function of sea clutter voltage is investigated 
through simulation. The results show that the number of breaking waves, which depends 
on the wind speed and the randomness factor of whitecaps, determine the shape of the 
probability density function of the clutter voltage. 

The randomness factor of the whitecaps causes a very 7 interesting change in the 
shape of probability density function. No familiar PDF can fit these forms. These 
probability density functions all meet the necessary 7 conditions that the breaking waves 
contribute to a tail in the high clutter voltage region. This explains the higher false 
alarm rate in detecting the presence of target return with a high resolution radar. 
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APPENDIX A SIMULATION PROGRAM OF BACK SCATTERING 



FROM WH I TEC APS 

define parameters and arrays 
REAL K 

REAL a: 8 SEED, S,RS( 200) 

DIMENSION U(20000),X(200),Y(200),C(200),AD(200) 
INTEGER D,NOS,NZ 

enter internal decorrelation time, seed for random generator, and 
whitecap randomness factor 
READ*, D,SEED,K 
define whitecap size and phase constant 
NZ= 128 
TH = 1 
Z= 1 

NN-2*NZ 

call uniform random number generator to get 2 random numbers 
as the seeds to be used in calling random generator for each 
zone size. 

CALL GGUBS(SEEDT83562.,NN,RS) 

PRINT*, NZ 

calculate scattering signal for each zone size. 

DO 1 1= 1,NZ 

determine number of zones to be updated. 

NS = ((I/D)+I)*2 

get phase and amplitude for zones to be updated. 

CALL GGUBS(RS(I)* 1006243., NS, U) 
update phase and amplitude. 

IC = (I/D)+ 1 
DO 2 J=1,IC 
JJ = J-1 
II = I-JJ*D 
X(II)=U(J) 

2 Y(II)= U(J + IC) 

A = 0. 
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C(I) = 0. 

calculate scattering signal. 

DO 3 L= 1,1 
IF(Z.EQ.O) THEN 

A = (- 1)**( L- 1 )*( 1 - K *X( L))*COS(TH) + A 
B = (- 1)**L*( 1 -K*X( L))*SIN(TH) + B 
ELSE 

A = (-l)**(L-l)*(l-K*X(L))*COS(3.14159 !l! K*(Y(L)-.5)) + A 
B = (-l)**L*(l-K*X(L))*SIN(3.14159*K*(Y(L)-.5)) + B 
END IF 

3 CONTINUE 

1 C(I) = A**2+ B**2 
E = 0. 

find the maximum scattering sample value. 

DO 4 I-l.NZ 
AD(I)= 1 

IF(C(I).GT.E) THEN 
E = C(I) 

END IF 

4 CONTINUE 
normalize all scattering signals 

DO 5 L= l.NZ 

5 C(L)= C(L) E 
plotting statements 

CALL RESET ('ALL') 

CALL TEK61S 
CALL PAGE (11. ,9) 

CALL HEIGHT (.2) 

CALL AREA2D(8.,6.5) 

CALL XNAMEfZONE NUMBERS',11) 

CALL YNAME('NORMALIZED RCSS',100) 

CALL HEADINfWHITE CAP SIMULATIONS', 100, 1.5,1) 
CALL GRAF(0, 16, 128,0, .1,1) 

CALL THKCRV (.01) 

CALL CURVE (AD,C,NZ,0) 
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CALL FRAME 
CALL HE1GHT(. 15) 

CALL XGRAXS(0.,L t 9.,8.,0,0,0.,0.) 
CALL YGRAXS(0.,1.,9.,S.,0,0,0.,0.) 
CALL BASALF('STAXD') 

CALL MIXALF('STAND') 

CALL REALXO(K,105,1.7,4.9) 
CALL RLMESSfK = S’, 3, 1.4, 5.5) 
CALL RLMESS('D = S’, 3, 1.4, 5. 2) 
CALL IXTXO(D,2,4.6) 

CALL ENDPL (0) 

CALL DOXEPL 

STOP 

EXD 
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APPENDIX B SIMULATION PROGRAM OF SEA SURFACE MODEL 



defining parameters and arrays 
REALMS SEED,DSEED,PSEED 

DIMENSION SE(9),PO(9 t 130),POIT(I30),NR(9),PSE(18),XP(9,I30) 
DIMENSION XTEP(130),YTEP(130),IP(9,10),CAP(9,10),IDX(9) 
DIMENSION YP(9,130) 
seed for random number generator 
SEED = 44567 

set up parameters to call poison random number generator 
TU= 100 
NUB =130 

get seeds to call poisson random number generator for nine areas 
CALL GGUBS(SEED*H 1 1000., 9, SE) 
get seeds to call uniform random number generator for locations 
CALL GGUBS(SEED*2981 1 1 17.,18,PSE) 
proceed breaking occurrence for nine areas 
DO 1 1 = 1,9 
DSEED = SE(I) 
get occurring times 

CALL POlSUB(DSEED,l.,0.,TU,NUB,POIT,N) 
two dimensional array to contain occurrence times 
DO 2 J = 1,N 
2 PO(I,J)= POIT(J) 

array to contain number of white caps in each area 
1 NR(I)= N 

loop to get positions for white caps in each area 
DO 5 1=1,9 
11 = 1 + 9 
KK = NR(I) 

get x and y position separately 
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CALL GGUBS(PSEED,KK,XTEP) 

PSEED = PSE(II) 

CALL GGUBS(PSEED,KK,YTEP) 
change references of position for sub-area 
FT = (I-.l) 

IF- FT/ 3 
IH = 0 

IF(I.EQ.2) THEN 
IH = 1 
END IF 

IF(I.EQ.5) THEN 
IH = 1 
END IF 

IF(I.EQ.8) THEN 
IH = I 
END IF 

I F( I.EQ.3) THEN 
IH = 2 
END IF 

IF(I.EQ.6) THEN 
IH = 2 
END1F 

IF(I.EQ.9) THEN 

IH = 2 

END IF 

DO 6 J = 1,KK 

XP(I,J) = XTEP(J)*5 + IH*5. 

6 YP(I,J)=YTEP(J)*5 + IF*5. 

5 CONTINUE 

loop to sample ocean surface in one resolution cell 
DO 2000 KKK= 1,50 
sample increment 0.2 sec 
SAMPT-KKK* 2+25. 
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scan all sub-areas 
DO 9 1-1,9 
1 1 F = 0 

check all white occurrence times in each sub-area 
DO 10 J= 1,NR(I) 

determine existing whitecaps and sizes at sampling time 
DIF= SAMPT-PO(I,J) 

1F(DIF.GT.0.) THEN 
1F(D1F.LE.(0.96)) THEN 
1 1 F = IIF+ 1 

store the size of each existing vvhitecap 
CAP(I,IIF) = DIF 
PRINT*, U,IIF,CAP(I,IIF) 

IP(I,IIF) = J 
END IF 
ENDIF 

10 CONTINUE 
IDX(I)= 1 1 F 
9 CONTINUE 
beginning of plotting 
CALL RESET ('ALL') 

CALL TEK61S 
CALL PAGE (11. ,11) 

CALL HEIGHT (.2) 

CALL AREA2D(8.,8.) 

CALL XNAME('X POSITION (M)S',14) 

CALL YNAME('Y POSITION (M)S',14) 

CALL H E AD I N( ' WH I TE CAP MODELS', 100, 1.5,1) 
CALL GRAF(0,5, 15,0,5, 15) 

DO 15 1=1,9 
15 PRINT*, I,IDX(I) 

loop to plot the sizes of existing white caps in it's location 
DO 11 1= 1,9 

IF(IDX(I).EQ.O) GO TO 11 
DO 13 J = 1,IDX(I) 
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XI = X P( I , I P( I , J ))/ 2. 

Y1 = YP(I,IP(I,J))/2. 

XN=CAP(U)/.0075 

CALL REA LN 0( C AP( I , J ) , 1 05 ,X 1 , Y 1 ) 

CALL REALNO(TU, 105,6,9.2) 

A = .9 

CALL REALNO(A, 105,6,8.9) 

13 CALL REALNO(SAMPT, 105, 6,8.6) 
11 CONTINUE 
CALL FRAME 
CALL THKCRV (.01) 

CALL ENDPL (1) 

2000 CONTINUE 
CALL DOXEPL 
STOP 
END 



subroutine to call poisson random number generator 

SUBROUTINE POISUB(DSEED,RLAMAX,RLAMIN,TU,NUB,P,N) 

REAL'S DSEED 
EXTERNAL A 
DIMENSION P(500) 

CALL GGNPP(DSEED !:s 234567,0.JU,NUB,A,RLAMAX,RLAMIN,0,N,P,IER) 

RETURN 

END 

FUNCTION A(TL) 

A = 1 

RETURN 

END 
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APPENDIX C SIMULATION PROGRAM OF TOTAL SCATTERING SIGNAL 
AND PROBABILITY DENSITY FUNCTION 



define arrays and parameters 
REAL-8 SD1,SD2,SEED 

DIMENSION SCM(20000),TFI(20000),PH(1500,128),CM(1500,128) 
DIMENSION XC( 1 50),I P(50),NSIZE(50),AD( 1 28),YSEQ( 1 500) 
DIMENSION TSEQ(0:1500) 

REAL K 
INTEGER D 

enter seeds for uniform generator and poison generator, roughness 
factor, and internal decorrelation time 
READ-, SD1,SD2,SEED,K,D 
XN2= K 

call subroutine to get a sequence of starting times, locations 
of breaking waves, and 128 samples of scattering for each whitecap 
over the whole observing time. 

CALL FINAL(NT,TSEQ,YSEQ,PH,CM,XC,K,D.SEED) 
call subroutine to get 20000 samples of scattering from unbroken 
water. 

CALL NOCAP(SDl,SD2,SCM,TH,AE,XE) 
sum the total scattering signal from unbroken water and breaking waves 
at 20000 sample times. 

DO 1 1= 1,20000 
set sample time. 

SAMPT= (.0075/2.) + (I- 1)*D*. 0075 
set initial existing breaking wave number at this time zero 
I IF= 0 

find whitecaps existing at this sample time 
DO 501 J = 1 ,NT 
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second 



DT = SAMPT-TSEQ(J) 

IF(DT.GT.O) THEN 
IF(DT.LE..96) THEN 
IIF= IIF+ 1 

determine zone size of whitecap 
1 1 S = DT/.0075 
NSIZE(IIF)= IIS 
store white cap sequence number 
I P( 1 1 F) = J 
END IF 
END IF 

501 CONTINUE 

get real part and imaginary part of scattering signal from 
unbroken water 

AX=SCM(I)*COS(TH(I)) 

AY=SC.M(I)*SIN(TH(I)) 

correct phase of scattering signal from breaking wave according to 
it's distance from the outer boundary of resolution cell, then add 
them to the scattering from unbroken water. 

DO 503 KK= 1,IIF 

THD = ( 15.-YSEQ(IP(KK)))*2*2*3.14159/.03 
PH(IP(KK),NSIZE(KK)) = PH(IP(KK),NSIZE(KK)) + THD 
AX = AX + C\l(IP(KK),NSIZE(KK))*COS(PH(IP(KK),NSIZE(KK))) 
503 AY = AY + CM(IP(KK),NSIZE(KK))*SIN(PH(IP(KK),NSIZE(KK))) 
get clutter voltage 

SCM(I)= SQRT(AX**2 + AY**2) 
calculate phase 

CALL PHASE(AX,AY,TH(I)) 

TEP= XN2 1 CONTINUE 
NN 1 = D 

DO 201 1 = 1,20000 
201 TH(I)= SCM(I)**2 
FID = 1 

calculate probability density function of sea clutter voltage 
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CALL PROB(20000,SCM,XT,XNT,XX2,FID) 

XX2= TEP 
FID = 0 

calculate probability density function of clutter cross section 
CALL PROB{ 20000, TH, XT, XXI, XX2, FID) 

CALL DOXEPL 

STOP 

EXD 



subroutine to get breaking starting times, locations, and scattering 
in the growing process of whitecap. this is the same as program in 
chapter 3 except calculation of scattering signal 

SUBROUTIXE FIXAL(XT,TSEQ,YSEQ,PH,CM,XC,K,D,SEED) 

REAL K 
INTEGER D 

REAL'S SEED,DSEED,PSEED,CSEED 

DI MEN’S I OX SE(9),PO(9,200),POIT(200),XR(9),PSE{ 1 S),XP(9,200 
DIMENSION’ XTEP(200),YTEP(200),1P(100),XS1ZE(100),CM( 1500, 128) 
D I M EXSIOX YSEQ( 1 500).PH( 1 500, 1 28),XC( 1 500) 

DIMENSION’ YP(9, 200), TSEQ(0: 1500) 

TU = 750 
XUB = 900 

CALL GGUBS(SEED*1 1 1 1000., 9, SE) 

CALL GGUBS(SEED*29S1 1 1 17..18.PSE) 

DO 1 1= 1,9 
DSEED = SE(I) 

CALL POISUB(DSEED,1.,0.,TU,XL'B,POIT,X) 

DO 2 J = l.X 
2 PO(I,J)= POIT(J) 

1 XR(I) = X 
DO 5 1-1,9 
11=1 + 9 
KK = XR(I) 

PSEED = PSE(I) 

CALL GGUBS(PSEED,KK,XTEP) 
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PSEED= PSE(II) 

CALL GGUBS(PSEED,KK,YTEP) 

FT= (I-.l) 

IF= FT/ 3 
IH = 0 

IF(I.EQ.2) THEN 
IH = 1 
END IF 

IF(I.EQ.5) THEN 
IH = 1 
END IF 

IF(I.EQ.S) THEN 
IH = 1 
END IF 

IF(I.EQ.3) THEN 
IH = 2 
END IF 

IF(I.EQ.6) THEN 
IH = 2 
END IF 

I F( I .EQ.9) THEN 

IH = 2 

END IF 

D06 J=1,KK 

XP( I ,J) = XTEP(J)*5 + IH *5. 

6 YP(I,J)= YTEP(J)*5 + IF*5. 

5 CONTINUE 

sort occurring times of all nine sub-areas in ascending order 
CALL TSORT(PO,NR,TSEQ,NT,YP,YSEQ) 

DO 601 1= l.NT 

set seed for random generator and call subroutine to calculate 
back scattering signal from breaking wave 
CSEED = TSEQ(I)*1 443567 
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CALL ZON"EM(D.CSEED,K,CM,PH,I) 

601 CONTINUE 

normalize scattering voltage from breaking by 1.224745 as the mean 



DO 605 1= 1,NT 
AVER = 0 
DO 606 J= 1,128 

606 AVER = CM(I,J) + AVER 
AVER = AVER/ 128. 

DO 607 K = 1,128 

607 CM(I,K)= CM(I,K)* 1.224745/AVER 
find maximum sample 

CALL M AX(CM ,1 ,XC) 

605 CONTINUE 
RETURN 
END 



subroutine to use poisson random number generator to get 
occurring times of breaking wave. 

SUBROUTINE POISUB(DSEED,RLAMAX,RLAMIN,TU\NUB,P,N) 

REAL* 8 DSEED 
EXTERNAL A 
DIMENSION P(500) 

CALL GGNPP(DSEED*234567,0.,TU,NUB,A,RLAMAX,RLAMIN,0,N\P.IER) 

RETURN- 

END 

FUNCTION A(TL) 

A = .100001 

RETURN- 

END 



subroutine to sort starting times of breakings from all sub- 
areas in ascending order. 
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SUBROUTINE TSORT(PO,XR,TSEQ,NT,YP,YSEQ) 
DIMENSION PO(9,200),NR(9),TSEQ(0:1500),D(1500),PT(1500) 
DIMENSION YT( 1 500), YP(9,200), YSEQ( 1 500) 

NT = 0 
DO 1 1=1,9 
1 NT = NT + NR(I) 

X = 1 

DO 9 1= 1,9 
DO 10 J=1,NR(I) 

YT(N)= YP(I,J) 

PT(N)= PO(I,J) 

10 N = X + 1 
9 CONTINUE 
TSEQ(0) = 0 
DO 8 L=1,XT 
DO 4 K-l.NT 
D(K) = PT(K)-TSEQ(L-1) 

IF(D(K).LE.0.) THEN 
D(K)= 750 
END IF 

4 CONTINUE 
TMIN= 750 
NP= 1 

DO 6 1= 1,NT 
I F(D(I ). LE.TM I N) THEN 
TMIN= D(I) 

NP= I 
END IF 

6 CONTINUE 
YSEQ(L) = YT(NP) 

8 TSEQ(L)= PT(NP) 

RETURN 

END 
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subroutine to calculate scattering signal from breaking wave 



SUBROUTINE ZOXEM(D.SEED,K,CM,PH,ID) 

REAL K 

REAL-8 SEED,S,RS(200) 

D I M ENS I ON U(20000),X(200),Y(200),CM ( 1 500, 1 28), PH( 1 500, 1 28) 
DIMENSION AD(128) 

INTEGER D,NOS,N’Z 
NZ= 128 
XN = 2*NZ 

get 2 seeds for each zone number 

CALL GGUBS(SEED- 100000000., N'N,RS) 

DO 1 1= 1,NZ 

compute number of random numbers needed for scattering signal at 
each zone number 
NS = ((I D)+l)*2 
get those random numbers 

CALL GGUBS(RS(I)* 100000000., NS, U) 
update phase and amplitude in zones according to internal decorre- 
lation time and current whitecap size. 

IC = (I D)+ 1 
DO 2 J= 1,IC 
JJ = J-1 
II = I-JJ-D 

IF(II.EQ.O) GO TO 7 
X(II)=U(J) 

2 Y(II)= U(J + IC) 

7 A = 0. 

B = 0. 

compute real and imaginary part of scattering from breaking wave 
DO 3 L= 1,1 

A = (-l)**(L-l)*(l-K*X(L))*COS(3.14159*K*(Y(L)-.5)) + ^ 

3 B = (-l)* ! r*(l-K*X(L))*SIN(3.14159*K !S: (Y(L)-.5)) + B 
compute phase 
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CALL PHASE(A,B,DC) 
PH(ID,I)= DC 
compute clutter voltage 
1 CM(ID,I) = SQRT(A**2 + B**2) 
RETURN 
END 



subroutine to calculate phase of a complex number 



SUBROUTINE PHASE(A,B,DC) 
DC = ATAN(B'A) 

IF(A.LE.O) THEN 
DC=DC + 3. 14159 
END IF 
RETURN 
END 



subroutine to Find the maximum number in an array 
SUBROUTINE MAX(CM,I,XC) 

DI MENS I ON CM( 1 500. 1 28),XC( 1 500) 

XC(I) = 0. 

DO 1 11= 1,128 
IF(XC(I).LE.CM(I,II)) THEN 
XC(I)=CM(I,II) 

END1F 

1 CONTINUE 
RETURN 
END 



subroutine to calculate scattering from unbroken water 



SUBROUTINE NOCAP(SDl,SD2,SCM,TH,AE,XE) 
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REAL'S SD1.SD2.PSEED 
DIMENSION SCM(20000),TH(20000) 
call Gaussian random generator to get random numbers for I and Q 
channel to get Rayleigh clutter voltage 
CALL GGNPM(SD1 *10023, 20000, SCM) 

CALL GGNPM(SD2*10255, 20000, TH) 

XE = 0. 

AE = 0. 

DO 1 1= 1,20000 
TCM = SCM(I) 

TTH = TH(I) 

SCM(I)= SQRT(TCM**2 + rTH**2) 

PSEED = I 

AE = AE + SCM(I) 

IF(XE.LE.SCM(I)> THEN 
XE= SCM(I) 

ENDIF 

1 CONTINUE 

call uniform generator to get random phase between plus minus n 
CALL GGUBS(PSEED, 20000, TH) 

AE = AE, 20000 
DO 12 1= 1,20000 
TH(I) = TH(I)*2*3. 14159 
12 SCM(I)= SCM(I)*.15;AE 

scale the clutter voltage by 0.15m as the mean for unbroken water. 
XE = XE ! M5/AE 
AE= .15 
RETURN 
END 



subroutine to calculate average in an array 



SUBROUTINE AVERA(AV,CC, II, NJ) 
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DIMENSION' CC( 1500, 128) 

c=o. 

DO 1 1= 1,NJ 
1 C = C + CC(II,I) 

AV = C/NJ*1. 

RETURN 

END 



subroutine to calculate the probability density function from 
20000 independent scattering samples 



SUBROUTINE PROB(SZ, SAMPLE, NT.NN1.XN2, FID) 
REAL K 
INTEGER SZ.D 

DIMENSION SAMPLE(20000),P(500),HX(500) 

REAL RA( 500), RA 1(500), RA2( 500) 

XS = 0. 

DO 100 1= 1,SZ 
IF(XS.LE.SAMPLE(I» THEN 
XS= SAMPLE(I) 

ENDIF 

100 CONTINUE 
XM = 0. 

fine the mean of samples 
DO 101 1= l.SZ 

101 XM = XM + SAMPLE(I) 

XM-XM/SZ 

set initial probability to be zero 
DO 105 1= 1,500 
105 P(I)=0. 

set the interval of clutter voltage to count number of samples 
located in them. 



IF(FID.EQ.O) THEN 
XD= .016 
ENDIF 



80 



JK = 0 

DO 102 1= 1,SZ 

determine which interval does the sample locate. 

SAP = SAMPLE(I), XD 
I SAP = SAP 
1 1 1 = 1SAP+ 1 

IF(III.GT.500) GO TO 102 
JK = JK+ 1 
P(III)=P(III)+1 

102 CONTINUE 
RSZ=SZ 

set initial value of total probability to be 0. 

TP = 0. 

set parameter for Rayleigh distribution according to sample mean 
AAR = XM '(SQRT(3. 14159/2.)) 
set initial value of total probability of Rayleigh distribution 
to be zero. 

TP1 = 0. 

compute probability density and total probability 
DO 103 1= 1,500 
HX(I)= I*XD-.5*XD 
P(I)= P(I)'(RSZ*XD) 

RA(I) = (HX(I).'AAR**2)*EXP(-1*HX(I)**2/(2*AAR**2)) 
TP1 = TP1 4- RA(I)*XD 

103 TP = TP + P(I)*XD 
PRC = JK/RSZ 

PRINT*, 'TP', TP, 'PERCENT', PRC 
statements to plot the probability density function 
CALL RESET ('ALL') 

CALL TEK618 
CALL PAGE (10.5,9.5) 

CALL HEIGHT (.15) 

CALL AREA2D(7.,7.) 

IF(FID.EQ.O) THEN 
CALL XNAME('RCSS',3) 
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CALL YNAMECPROBABLITY DENSITY (1/M**2) S', 28) 

ELSE 

CALL XXAMEf V S', 3) 

CALL YNAME('PROBABLITY DENSITY (1/M)S',24) 

END IF 

CALL HEADIN('PROBABILITY DENSITY FUNCTIONS', 28, 1.5,1) 
IF(FID.EQ.O) THEN 
CALL GRAF(0,.6,6,0,3.6,36) 

ELSE 

CALL GRAF(0,.8,8,0,.24,2.4) 

ENDIF 

CALL THKCRV (.01) 

CALL MARKER(15) 

CALL SCLPIC(.IS) 

CALL CURVE (HX,P, 500,-1) 

CALL HEIGHT(.15) 

CALL XGRAXS(0.,1.,9.,8.,0,0,0.,0.) 

CALL YGRAXS(0.,1.,9.,8.,0,0,0.,0.) 

CALL BASALF( STAND') 

CALL MIXALF('STAND') 

IF(FID.EQ.O) THEN 

CALL RLMESSfMAX RCSS',7,6.3,6.43) 

PRINT*, XS,'Q' 

CALL REALNO(XS, 105,6.7,5.73) 

CALL RLMESS('MEAN RCSS',8,6.15,6.1) 

CALL REALNO(XM, 105,6.7,5.43) 

ELSE 

CALL RLMESS('MAX V S', 7, 6.3, 6.43) 

CALL REALNO(XS, 105,6.7,5.73) 

CALL RLMESS('MEAN V S', 8, 6.15, 6.1) 

CALL REALNO(XM, 105, 6.7, 5.43) 

ENDIF 

CALL RLMESS('TOTAL P S', 8, 6.15, 5.77) 

CALL REALNO(TP, 105, 6.7, 5.13) 
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CALL RLMESS(TOTAL WHITE CAPSS',16,5,7.5) 
CALL INTNO(NT, 6.8, 6.65) 

CALL RLMESS('DSM, 7.1, 7.15) 

CALL INTNO(NNl, 6.8, 6.35) 

CALL RLMESS('KS', 1,7.1, 6.8) 

CALL REALN0(XX2, 105, 6.7, 6.05) 

PRINT*, TP1 
CALL EXDPL (0) 

RETURN 

EXD 
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